www.gusucode.com > Weighted Differential Evolution Algorithm (WDE) > Weighted Differential Evolution Algorithm (WDE)/cc_wde/PointCloudProcessing/my_tin_simplifier.m
function [out,fnc,v]=my_tin_simplifier(X,mydata) n0=mydata.n; xtest=mydata.xs; ytest=mydata.ys; ztest=mydata.zs; fnctest=mydata.fnc; N=size(X,1); out=rand(N,1); % memory for i=1:N sol=X(i,:); sol=reshape(sol,n0,2); % planimetric locations of tie points zsol=fnctest(sol(:,1),sol(:,2)); v=[sol zsol]; fnc = scatteredInterpolant(v(:,1),v(:,2),v(:,3),'natural'); z0=fnc(xtest,ytest); err=sum(abs(z0-ztest)); out(i)=err; end